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Flavour singlet contributions to the nucleon spin are elusive due to the fact that they 
cannot be determined directly in experiment but require extrapolations to the small x 
region. Direct calculations of these contributions are possible using Lattice QCD, however, 
they pose a significant computational challenge due to the presence of disconnected quark 
line diagrams. We report on recent progress in determining these sea quark contributions 
on the lattice. 



1 Introduction and Results 

The distribution of the spin of the proton among its constituents has long been a topic of 
interest. The total spin can be decomposed into the contribution from the quark spins, AS, 
the quark orbital angular momenta, L^, and the gluon total angular momentum Jg [T], 

I = ^AS + L^ + J<„ (1) 

where AS = Au + Ad + As (heavier quarks are normally neglected). In this work, Aq (q = 
u, d, s) denotes the combined spin contribution of the quark and the antiquark. Using Lattice 
QCD, one can determine the Aq from first principles through the axial-vector matrix element, 

— (7V,s|g7^75jg|7V,s) = ^ (2) 

where m^r is the mass of the nucleon with spin (s^ ~ ^1)- Thus, one can construct the axial 
charges, as = Au-Ad = gA, as = Au + Ad-2As and ao(Q^) = Au + Ad+As = AT,{Q'^). AS 
acquires a scale dependence, Q^, due to the axial anomaly. The axial- vector matrix element 
is related to the first moment of the quark helicity distributions. The second moment, (a;)A, 
and the second moment of the transverse helicity distribution, {x)s, can also be calculated on 
the lattice, higher moments are more challenging. Furthermore, the total angular momentum 
of quark q, Jq = ^Aq + Lq {J2q Jq = + -^vO '^^'^ obtained from the generalised form 

factors, v42o(Q^) and i?|Q((5^), which parameterise the matrix element of the energy-momentum 
tensor for momentum transfer, Q. For a review of recent Lattice results of Jq, (a::)A and {x)s 
see [2 13] and references therein. 

In these proceedings we focus on Aq and As in particular. Lattice results for As have an 
important role to play in constraining fits of polarised parton distribution functions (PDF). 
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The spin structure function of the proton and neutron, Q^), is measured in deep in- 

elastic experiments. The first moment is related to the axial charges via the operator product 
expansion. To leading twist: 



where Cs and Cns are the singlet and non-singlet Wilson coefficients, respectively. Model 
assumptions are made in order to extrapolate 17"'^ (a;, Q^) from the minimum x = 0.02 accessible 
in experiment down to a; = 0. 03 is known from neutron /3-decay, while, assuming SU{3) flavour 
symmetry, og can be obtained from hyperon /3-decays. Thus, in combination with r^'"((5^), 
oo and the Aq can be deduced. For example, HERMES find As(5GeV^) = 5(00 — as) — 
— 0.085(13)(8)(9) [1]. However, if the range of x in the integral in Eq. [s] is restricted to the 
experimental range, x > 0.02, As is consistent with zero, indicating the large negative value 
arises from model assumptions in the low x region. 

Semi-inclusive deep inelastic scattering (SIDIS) experiments offer a direct measurement of 
the Aq using pion and kaon beams. Results from COMPASS show the strangeness contri- 
bution is consistent with zero down to x = 0.004 [S]. A naive extrapolation to x = gives 
As = —0.02(2) (2), while using the parameterisation of De Florian et al. (DSSV) gives 
As = — 0.10(2)(2). Present measurements via SIDIS are limited by the knowledge of the quark 
fragmentation functions, to which As is particularly sensitive. Another possibility of directly 
determining As combines z/p. Dp and parity- violating ep elastic scattering data [7 . Here, the 
MicroBooNE experiment will enable errors to be significantly reduced. 

Considering the lattice approach, simulations are performed at finite volume (V^) and lattice 
spacing (a) and typically with u and d quarks with unphysically heavy masses (mg). Physical 
results are recovered in the continuum (a — 0) and infinite volume {V — > 00) limits at physical 
quark masses. The u, d and s quark masses used in simulations are normally expressed in 
terms of the pseudoscalar meson masses they correspond to {Mpg oc iriq). Developments in 
algorithmic techniques and computing power mean typical simulations now involve lattices with 
a - 0.05 - 0.1 fm, F - 2.5 - 6 fm and Mps{u/d) - 200 MeV. 

For any lattice prediction, the size of the main systematic errors, discretisation effects, fi- 
nite volume and so on, must be investigated thoroughly. However, the systematic uncertainties 
should always be compared to the inherent statistical error, where for some quantities the latter 
dominates. As is one such example. Evaluating {N\s"f^j5s\N) involves calculating a discon- 
nected quark line diagram (for the strange quark). These types of diagrams are computationally 
expensive to calculate as they involve the quark propagator from all space-time lattice points 
to all points. In the past these diagrams were often not calculated and differences of quantities 
were quoted, for which the disconnected contribution cancelled assuming isospin symmetry, for 
example, qa = Au — Ad, Ju-d and {x)Au~Ad- However, methods have been developed which 
enable the disconnected contributions to be calculated [8], in particular, also As. Note that 
for Am and Ad a "connected" quark line diagram must also be evaluated in addition to the 
disconnected one. 

In the following we present results for the quark spin contributions to the proton generated 
on configurations with two degenerate flavours of sea quark^ with a ^ 0.072 fm and u/d quark 
masses given by M ps = 285 MeV. The quark action employed has leading order discretisation 

^The number of sea quarks refers to the number of flavours of quark fields included in the Monte Carlo 
generation of an ensemble of representative quark and gluon field configurations. 
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effects of O(a^), which are not expected to be significant for this value of the lattice spacing. 
Two lattice volumes, with spatial dimensions of 2.3 fm and 2.9 fm, were used and no significant 
finite volume effects were found. Although the u and d quark masses are unphysically heavy, 
we varied the quark masses in the range corresponding to Mps — 285 — 720 MeV and found 
no significant change in the results. This suggests our results may also apply to physical quark 
masses. Further details can be found in [9]. 
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^ lat 
^ycon 






DSSV XjYiiji 


DSSV 


u 


1.071(15) 


-0.049(17) 


0.787(18)(2) 


0.793 


0.814 


d 


-0.369( 9) 


-0.049(17) 


-0.319(15)(1) 


-0.416 


-0.458 


s 





-0.027(12) 


-0.020(10)(1) 


-0.012 


-0.114 


as 


1.439(17) 





1.105(13)(2) 


1.21 


1.272 




0.702(18) 


-0.044(19) 


0.507(20)(1) 


0.401 


0.583 


E 


0.702(18) 


-0.124(44) 


0.448(37) (2) 


0.366 


0.242 



Table 1: The bare (lattice) connected, Aq'^'^, and disconnected, Agjg*, quark spin contributions 
of the proton. Renormalisation to the MS scheme is performed for Q ~ \/YA GeV, to obtain 
Ag^^((5). The first error is statistical, the second is from the renormalisation. A comparison 
is made with the results of DSSV global fits to DIS and SIDIS data [TU] (see this reference for 
the uncertainties on the values) for the cases where x is restricted to [Xmin = 0.001 — >■ 1] and 
where x = [0 ^ 1] . 

Table [l] displays the results for As and the axial charges and compares them to results from 
a DSSV global analysis [TU]. We obtain a small, negative value for As, consistent with the 
result obtained from a truncated DSSV fit to experimental data. Notably, 03 is significantly 
lower than the experimental value of = gA = 1.2695(29). This is a general feature of lattice 
calculations of gA- A summary plot of recent simulations taken from is shown in Fig. [l] 
Over the range of u/ d quark masses available the results obtained using different lattice quark 
actions, volumes and lattice spacings, are fairly constant and lie approximately 10% below 
the experimental result. Simulations at smaller u/d quark masses are needed. Problems with 
excited state contaminations have also been suggested as a possible source of the discrepancy, 
see, for example, Ref. [11 . On the lattice an operator with the quantum numbers of the proton, 
creates all states of = ^ and one needs to ensure that the matrix element for the ground 
state has been extracted. Nonetheless, the uncertainty in the lattice determination of gA seems 
to be multiplicative, as shown on the right in Fig. [I] the ratio of gA/ f-K tends to the experimental 
value. 

Considering the underestimate of gA, we add a 20% uncertainty to our results and obtain a 
final value of 

As = -0.020(10)(4). 

The first error is statistical, the second is due to the systematic uncertainty, which dominates. 
Other groups using similar methods, albeit at heavier quark masses and without renormalisa- 
tion, obtain consistent results, see |12| and |13| . This is in contradiction to earlier exploratory 
work by, for example, the Kentucky group |14j . Note that this group have recently also calcu- 
lated flavour singlet contributions to Jq [T5] . 
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Figure 1: (left) A summary plot of — for different lattice simulations as a function of the 
u/d quark mass, expressed in terms of Mpg = rri^ [3]. (right) The ratio of gA/ f-n as a function 
of for different lattice spacings TB]. /3 = 5.2 — 5.4 corresponds to a ^ 0.083 — 0.060 fm. 
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